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ABSTR/rr 


We show how relatively simple, non-imaging, microwave radars can be used to measure ocean 
wave directional spectra from satellites. In the non-imaging radar measurement, die directional dis- 
tribution of tile ocean waves is obtained by scanning the radar in azimuth; the resolution of Fourier 
component waves in any given azimuth direction is a result of a Bragg type of wave front matching 
condition. The spectral distribution in wave number can in principle be measured using either short 
pulse (SP) or dual frequency (DF) techniques. 

In this paper, we analyze the SP and DF techniques in the frequency domain, in which is evi- 
dent a common dependence on a generalized fourth-order statistical moment of the surface scatter- 
ing transfer function. A solution for the fourth-order moment is obtained using physical optics in 
the high frequency limit appropriate to near-vertical, specular backscatter. The solution is expressed 
as an asymptotic expansion valid to the second-order in Gaussian wave statistics and the first order 
in non-Gaussian wave statistics. The moment solution gives, to the first order in wave steepness, an 
intrinsic electromagnetic modulation spectrum that is proportional to the large wave directional 
slope spectrum evaluated in the azimuth of radar look. The harmonic distortion inherent to the 
specular scatter measurement of the sea slope spectrum is found to be small in the incidence angle 
range of 8 to IS degrees. 

Detailed consideration is given to the measurement signal-to-noise problems specific to SP and 
DF techniques. It is shown that with suitable pulse integration, typical satellite measurement signal- 
to-noise ratios of ca. 0 dB result for the narrowband DF technique and -<-20 dB for the SP technique. 
Thus, while satellite measurement feasibility with either technique is indicated, the SP technique 
possesses a distinct advantage in terms of measurement signal-to-noise and contrast ratios. 
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1. Introduction 

'Traveling by air over the ocean, have you ever noticed how dearly visible were the patterns of 
the ocean waves when direct sunlight illuminated the surface? lask the question because I want to 
nuke the point that if these waves are visible to the naked eye, then study they ought tobe meas- 
urable by some remote sensing technique. In particular let us inquire if they cannot be measured by 
microwave radars which, as remote sensors, possess the distinct advantages of (a) providing their own 
source of illumination, (b) seeing through clouds and light rain and (c) being practically insensitive 
to the ionosphere. We will show that indeed microwave radars can measure the waves; moreover we 
show that the radars do not have to image the sea surface, if only statistical properties of the wave 
field are desired. The most useful statistical description of a random, homogeneous wave field is 
given by the directional wave height spectrum. In terms of the spatial properties of the wave field 
that are most easily observable from remote platforms, this is the two-dimensional, vector wave num- 
ber spectrum of the two-dimensional height field. 

Modem, numerical wave forecasting and analysis is based on wave spectral growth and propaga- 
tion models that require as inputs initial sea state (i.e., directional wave spectra) and future wind 
field specifications. Owing to inevitable model and specification errors, the wave forecasts will be 
in error. Presently, there are no viable means, either in-si tu or remote, of measuring directional ccean 
wave spectra in any kind of routine or global sense. If routine remote observations were to become 
available, we should be in a much better position to verify the wave forecasts currently being pro- 
duced and to improve upon them, through both model improvement and improvement of the initial 
state specification. 

We shall be considering observations of the sea surface in the case where the radar is pointing 
downward, away from the vertical, but not more than 1 5 degrees or so from the vertical. As is well 
known (Barrick, 1968a, 1974; Valenzuela, 1978), near-vertical microwave backscatter from the sea 
occurs by means of quasi-specular retroreflections from surface wave facets that are oriented normal 
to the radar line-of-sight. The backscatter can be described to various levels of approximation by 


3 


gefflhefriCBfoptics, physical optics, or ttigher-order high-frequency approximations (Jackson, 1974a); 
or alternatively, by composite surface theory (Brown, 1978). Basically though, the situation is one 
of specular reflection, and the radar backscatter problem can be compared to that of sun glint ob* 
servation in the case where the sun is anti-parallel to flut observer’s line ofsight. 

I am concentrating here on small angle, specular backscatter for two teaseiHi One is that Ibe* 
Ueve that the phfW K mechanisms involved in tire reflectivity modulation by tiler large waves are 
simpler than in the case of large-angle Bragg diffraction backscatter. At large angles there is, in addh 
tion to a purely geometrical tilting effect, a very strong sensitivity of the electromagnetic (Cm) mod- 
ulation to the details of the hydrodynamic modulation of the centimetric Bragg-diffracting wavelets 
by the atmospheric and large ocean wave flow fields, While the hydrodynamic contribution to the 
large-angle reflectivity modulation has been treated theoretically (Alpers and Hpselmann, 1978; 
Valenzuela and Wright, 1 979), theory is hard pressed to deal with the hydrodynamic complexities 
and actual field observations (Valenzuela and Wright, 1979; Wright et aL, 1980). In the specular 
scatter regime, however, there is no special sensitivity to a particular small water wavelength. Rithtr, 
tee entire ensemble of water waves is contributing to tee population of specularly reflecting facet* 
(excluding of course those wavelets smaller than the diffraction limit). For this large ensemble of 
waves, an assumption of Gaussian surface statistics seems to be a reasonable one, at least in the first 
order of approximation. 

The second reason for concentrating on the near-vertical is a more practical one. It is desirable 
to make the nadir angle as small as possible in order that the radius of the azimuth scan circle on the 
ocean surface not exceed the distance over which the wave field can be expected to be homogeneous 
(cf. Fig. 1 ). 

We have said that the radars need Hot actually image the surface to measure wave spectra. The 
idea behind the non-imaging radar measurement is simple: we let the broad cm phase front on the 
surface (i.e., broad compared to the correlation scale of the waves) function to isolate or resolve 
plane Fourier surface contrast waves travelling in (or contrary to) the direction of radar look. The 


4 


principle is similar to that of Bragg scatter, which holds that backscatter occurs only from those water 
waves having wave vectors "matched” to the em propagation vector. Specifically, the water wave 
propagation vector K must lie in the plane of incidence and must have the magnitude 

| K| -2*sin0 


where k is the em wave number and 0 is the angle of incidence. To illustrate the principle we con- 
sider a simple tilting model of reflectivity variation in which die variation is directly proportional to 
the large wave slope component in the plane of incidence. See Fig. 2. Then if Vf is the large wave 
dope, die radar modulation spectrum will be proportional to the large wave directional slope spec- 
trum evaluated in the azimuth of radar look, i.e., 


WK> « 



■Vjtg) e" 1 ^ * S dx 


|K| 2 F(K) 


where brackets < . . . > denote ensemble average, p is the unit horizontal vector, K/K, and F(K) is the 
spectrum of wave height 

The actual detection of the range-travelling surface contrast waves can be accomplished using 
either short pulse (SP) or dual frequency (DF) techniques. In the SP technique, very short pulses 
are used to resolve the surface wave structure in range. The fast-time or range modulation of the 
backsca tiered pulse caused by the spatial reflectivity variation can be detected by spectrum analysis 
of the backscattered signal envelope. (See Fig. 3.) The SP approach to possible satellite ocean wave 
measurement was apparently first suggested in a note by Tomiyasu (1971). Apart from the work of 
Jackson (1974b) and the present work, the SP approach has received no further consideration in the 
literature. 

The DF approach, first proposed by Ruck etaL (1972) (and cf. also Ransone and Wright, 1972; 
Hasselmann, 1972) has received considerably more attention (Plant, 1977; Alpers and Hasselmann, 
1978; Schuler, 1978; Plant and Schuler, 1 980). The DF technique was originally proposed as - 
and is still generally understood to be - a narrowband technique. Quasi-monochromatic waveforms 
that entirely fill the antenna beam spot on the surface are transmitted at two closely separated 
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microwave frequencies. An harmonic surface contrast modulation is detected through die interaction 
of the em “beat wave" with the surface. Jackson (1974b) showed that the narrowband DF technique 
k has a poor measurement SNR (signal-to-noise ratio) compared to the SP technique. The reason for 

the low SNR is that the contrast modulation spectrum is relatively broadband, having a bandwidth 

* i 

comparable to the bandwidth of the sea-slope spectrum (a bandwidth of the order of SO- 100 percent 

i ' 

of the peak frequency). 'Since the modulation spectrum P mo( j(K) is being detected only in a narrow H 

1 ■ J 

spectral band of width AK. ~ 2a/ L p where L p is the range extent of illumination, the detected signal 
(the covariance of powers at the two frequencies) will be small compared to the background noise j 

level (the variance of backscattered power at either frequency), i.e., the SNR ■ j 

Jackson ( 1 974b) showed that this low SNR problem inherent to a narrowband DF measurement could J 

be solved by using wideband signals - e.g., by using signals with bandwidths comparable to the sea- ij 

spectrum bandwidth. An illustration of a wideband DF radar system is given in Fig. 3. We note that -j 

A 

'A 

the receiving system represented in Fig. 3b, consisting of a two channel receiver and cross-correlator, 3 

was described some years ago by Parzen and Shiren (19S6). Parzen and Shiren considered both the 1 

SP and DF receiving systems from the viewpoint of detection theory. Independently of any particu- I 

J 

I lar transmitter, the SP and DF receiving systems were regarded simply as alternative systems for de- j 

tecting a modulation on a noisy carrier signal. In their analysis, Parzen and Shiren assumed that the I 

l input to both detection systems was a weakly modulated Gaussian noise process. This is actually a | 

1 realistic model of the sea-return process we are considering, and indeed a performance comparison 

of SP and DF systems can be made on the basis of such a model process (to be more specific, the 
surface impulse response can be modelled as a weakly modulated complex Gaussian process). 

Alpcrs and Hasselmann (1978) (hereafter referred to as A-H) performed an extensive analysis 
of the DF technique in its original narrowband (continuous wave) context. Their analysis dealt with 
large-angle backscatter and included hydrodynamic modulation. Rather than using wideband signals 
to improve the DF SNR as suggested by Jackson ( 1 974b), A-H instead introduced the idea of slow- 
time filtering the DF signal D, : (t) = H,(t)E^(t) where E, and F, are the backscattered fields at the 
two frequencies. The idea is that the modulation signal component of Dj,(t) will appear as a strong 
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Une in the comparatively broad Doppler spectrum of D )2 (t). A-H showed that by filtering D 12 (t) at 
the appropriate platform-motion-induced Doppler frequency, satellite SNRs typically of 0 dB could 
be obtained. 

In this work, we apply the A-H filtering scheme to the SP technique as well as to the DF tech- 
nique. In the SP technique, the slow-time filtering can be realized simply by integrating pulse re- 
turns in range bins (fast-time bins) that are adjusted so as to compensate for the apparent motion of 
the target ocean scene due to the motion of the platform. A schematic illustration of this kind of 
processing in a SP system is given in Fig. 4. (We remark that the integration, or filtering, is done to 
reduce the level of the random signal fluctuation arising from waveform coherency effects (Rayleigh 
fading). 

Following, we present a frequency-domain analysis of the SP and DF systems illustrated in Fig. 
3. This introduces the generalized fourth-ordei moment of the surface scattering transfer function. 
.Using the scalar physical optics integral solution for the transfer function, we obtain, first, a time- 
independent, far-zone solution for the moment in the high-frequency limit appropriate to near- 
vertical backscatter (geometrical optics limit). We proceed then from the far-zone solution to the 
more complicated time-dependent Fresnel zone moment solution that is required for an analysis of 
the systems employing the A-H slow-time filtering. We conclude with an SNR analysis of SP and 
DF techniques which indicates the feasibility of measuring ocean wave directional spectra from low 
earth-orbiting satellites. The analysis shows that the processing gain achieved with the A-H filtering 
scheme is approximately the ratio of the microwave carrier fiequency to the ocean-wave induced 
modulation frequency of interest. For X or Ku band radars, this ratio is on the order of +30 to +40 
dB. With this processing gain, we arrive at typical satellite measurements SNRs of ca, -10 dB for the 
narrowband DF approach and ca. +20 dB for the SP approach. 

2. Short-Pulse (SP) and Dual-Frequency (DF) Techniques 

Figure 1 depicts a pencil beam of radiation incident on the sea surface. The antenna beam axis 
makes an angle 6 to the vertical (z axis) and lies at an azimuth 4> relative to some (x t , x 2 ) reference 
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axes in the horizontal plane. The unit vector £ * (cos4>, sin 4) is defined as pointing along the azimuth 
ray, or in the direction of increasing surface range. The slant range from the radar to points on the sur- 
face z ■ f(x), x ■ (x j , x 2 ), is denoted r, and the range to the center of the beam spot is denoted r 0 . 

We shall be considering antenna beam widths typically of a few degrees and illuminated areas Aq 
typically of several kmin lateral and range extents. 

We consider the backscatter of any finite duration but otherwise arbitrary transmitted wavefomu 
f(r). Here we use r to denote “fast” signal.time, i.e., that time having a scale commensurate with the. 
pulse duration, or reciprocal pulse bandwidth (typical scale of microseconds). Fast time r is to be 
distinguished from “slow” time t, which scales with the interpulse period or reciprocal Doppler band- 
width (typical scale of milliseconds). The incident field at the surface is given by the spherical wave 
e^vr, x) = G w (x)f(r-r/c)/r, where G is the antenna power pattern on the surface. The Fourier 
transform (FT) of the incident field is 

E lnc (k,x) = G % (x)e- ikr E 0 (k) (la) 

where E 0 (k) is the waveform FT at the surface, 

E 0 (k) « (c/2 w)J*c -ikCr l f(r)/r 0 1 dr (lb) 

where k 8 v/c is the wave number, v is the frequency and c the speed of light; and where in the de- 
nominator we have set r 8 r 0 since, for the few degree beamwidths and large ranges we are consider- 
ing, r varies only very little over the beam spot. The backscattered field e ( (r,t) is a duration-limited 
random process in fast time r (the duration being determined by the round-trip travel time T r of an 
impulse across the beam spot on the surface) and is assumed to be a stationary random process in 
slow time t. Stationarity in t follows if wo«ssume that the radar is translating at a uniform velocity 
over a homogeneous sea. For simplicity, we will ignore the effect of discrete pulsing and sampling 
and treat e s (r,t) as a continuous process in t. Let the FT of e, with respect to r be denoted as 
E,(k, t). Then E, is related to the incident wave form FT E 0 (k) as 


E,(k,t) 8 S(k,t)E 0 (k) 


( 2 ) 
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where S(k,t) is the surface scattering transfer function for backscatter at a frequency v ■ kc. In 
principle, S(k,t) is tire solution to the complex boundary value problem for a unit amplitude har- 
monic incident wave in the case where the scattering geometry is frozen at some time t Note that 
we are treating S as a scalar quantity; this is permissible since them is only a very weak polarization 
dependence in the specular backscatter regime (cf. Jackson, 1974a). Also note that the look angle 
0,4> dependence is implicit in S. 

A frequency-domain analysis of the SP and DF systems diagrami'***! in Fig. 3 is given in Appen- 
dix A. The SP and DF systems of Fig. 3 will, for convenience of identification, be referred to as 
“modulation noncoherent” systems, as they are insensitive to the phase of die signal modulation. 

In contradistinction, those systems which are sensitive to the phase of the modulation, such as the 
system of Figure 4, will be referred to as ‘‘modulation coherent." Since die idealized SP and DF 
systems double detect using square law detectors, their outputs will depend on fourth-order products 
of the surface transfer function S(k,t). From Eqs. A2 and A8, it is seen that the ensemble average 
outputs of the SP and DF systems can be expressed in terms of different operation 1 , on the same 
“generalized” fourth-order statistical moment of S(k,t), viz., 

M(k,k',x) - <S(k)S*(k-x)S*(k')S(k'-K)> (3) 

where < . . . ) denotes ensemble average and * complex conjugation. Note that in (3) we have omitted 
the explicit time dependence in S(k,t). If die difference wave number A k is defined as 

Ak*k-k' (3*) 

then M can also be expressed as a function of Ak, M ■ M(k, Ak,*). In the following, we will be work- 
ing under the assumption that the pulse bandwidthsfor discrete difference frequencies) are not exces- 
sive; hence we require that 

«c/k« ) and Ak/k « I (3") 

In practice (3") will always be well satisfied; e^.. a 1 percent bandwidth would represent quite a large 


bandwidth. 


3. Specular Scatter Solution for the Moment M(k, Ak,<t) 

The physical optics (PO) integral solution for the transfer function can be written as (cf. Beck* 
ir.ann and Spizzichino, 1963; Weissman, 1973) 

S(M) " 7^J G fc-¥*>®^H2ki< s ,t)ida w 

wltere G is the illumination pattern; V is the platform velocity; j ■ (Xj , x 2 ) is the heiizontal coordi- 
nate vector of surface points in a surface-fixed reference frame; dx ■ dx ( dx 2 is the element of hori- 
zontal area; and where the integration extends over the illuminated area. The origin of the coordinate 
system x ■ Q is taken to coincide with the center of the beam spot at t * 0. For small t, the Fresnels 
approximation to the phase is valid, 

-2kr(x,t) * -2kr 0 ♦ 2kcos0f(x) - 2ksin0 g • (x - Vt) 

+ (k/r 0 ){| x- Vt| 2 -sin 2 0l£ *(x- Vt)l 2 + f 2 (x)} (5) 

where again f(x) is the random wave height and £ is the unit vector in the horizontal pointing in the 
direction of increasing range. In the analysis of the modulation noncoherent systems of Fig. 3, we 
find that a far zone approximation to the phase is adequate; the quadratic Fresnel term in (S) can be 
dispensed with without sacrificing any of the essential physics of the problem. Thus, neglecting the 
Fresnel term (and ignoring the time dependence which is now trivial without the Fresnel term), we 
have for the far zone approximation to the PO integral: 


Sik) * 


ilte^ 2kT o f 

— — lG(x)exp{i2klcos0f(x)- sin0 g • xj J dx 


( 6 ) 


We calculate the fourth moment, Eq. (3), in the geometrical optics (GO) limit appropriate to near 
vertical, specular backscatter. The GO limit of the moment (3) can be calculated in a number of 
different ways according to various methods that can be found in the literature dealing with the 
second moment <|S(k)| 2 > (average backscattered power). Of these methods, Barrick's ( 1 968b) is 
the most direct. Other methods that start with the PO integral include those of Beckmann and 
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Spizzichino (1963) and Kodii (1966: cf. Weinman, 1973). Of course, it b possible to dispense with 
PO altogether and start ab initio with GO, e.g., using methods such as Lynch and Wagner's (1970). 
But, as we still have to account for phase, there is rea’’y no advantage to be gained in using GO. Pro- 
ceeding according to Barrick’s( 1968b) method, we substitute the PO integral (6) into (3). Writing 
the product as a four-fold integral, and interchanging onseir ble average and integration operations, 
we have: 

M(k,Ak,x) - J" • • jGfx 1 ) . . . Gfk 4 ) 

• ^exp{i2kcosdlk(f 1 - f 2 - f 3 + i 4 ) ♦ Aktf 3 - f 4 ) 

♦Ktf 2 -! 4 ))}) • exp {-2sin0 £ • (k(x* - x 2 

-x 3 + x 4 ) + Ak(x 3 - 3 4 ) + k(s 2 - 3 4 )1| dx 1 . . . dj 4 (7) 

where we have approximated the product of the wave numbers k(k-x) ... by k 4 (which b permis- 
sible in view of (3*)), and where we have let f* ■ f(x l ), etc. For large k (i.e., for large *ms phase 
variation 2kcos0<f 2 >**), it b apparent that appreciable contributions to the moment are made only 
in the neighborhoods of two sets of stationary points, viz., 

S, : x 1 -x 2 , x 3 “x 4 

S 2 : x 1 "* x 3 , x 2 “x 4 (8) 

Except where x 1 ■ x 2 ■ x 3 ■ x 4 , the sets 5, and S 2 are dbtinct and yield distinct contributions 
to the moment, i.e., M ■ M, ♦ M 2 , where and M 2 respectively derive from integrations over the 
small volumes surrounding 5 } and 5 2 . Consider the integration over the volume containing 5 j . let 

u“ 2kcos0(x l -x 2 ) 
v ■ 2kcos0(x 3 - x 4 ) 
w ■ x 2 - x 4 


(9a) 


Now expand the height difference* in the neighborhoods of the stationary points S x in a Taylor 
series. 


It 


2kcos0(f 1 - f 2 ) ■ vf l *H + 0(u 2 /k) 

2kcosfl(f 3 - f 4 ) ■ v ♦ Of^/k) (9b) 

where V ■ (d/xj, d/dx 2 ). Then the contribution M t to M for large k becomes 


Mj ~ {ITqCos^B )" 4 J A(w) exp(-i2ksin0 £-w)dw 
X ^x)" 4 1 1 {exp{i|Vf 2 *si - vf 4 *^ 

♦ 2kcos 0(f 2 - f 4 ) + 0((u 2 + v 2 )/k) + OCAkv/k)]} ) 

• exp j i l land g-(u - v) 0(Akv/k» } du dv J 


( 10 ) 


where in (10) we have already taken the large k limit in the products of the gains, i.e., G(x' ) . . . Gfx 4 ) 

I ~ G 2 (x)G 2 (x + w), and where we h;ive integrated over the remaining space variable x, letting A(w) 

I ' “ 

represent the convolved two-wey gain pattern. 


i 

A(w) ■ 



G 2 (x)G 2 (x + w)dx 


(ID 


(Since f(x) is a homogeneous random process in x, the average ( • • • > is independent of **ie absolute 
position vector x ■ x 4 , so thu the x integration only effects the gain function.) 

Now consider the limit of the factor in braces as k -* «. Lct^ denote the random six-vector, 

J-W 2 .?* 4 .* 2 .* 4 ) (12a> 


and let t denote the associated six-dimensional characteristic vector, 


t » (jj. -v, 2*cos0, -2kcos 0) (12b) 

Then, the characteristic function of £ is defined as 

(t;w)» (e 1 !’*^- Je*i’i P t (|;w)d| (13) 

where p { is the pdf (probability density function) of Since f(x) (and hence |) is a stationary proc- 
ess, ti e pdf and characteristic function depend only on the separation, or lag vector w * y - x 4 . 



i 

i 

i 
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From the definitions (1 2, 13) we have for M t in the high-frequency, GO limit, at least formally, 

lim M j s M| (k) = (2TQCOS0)" 4 I A(w) exp(-i2xsin0 g*w)dw 
k-*» J 

X (ttJr)" 4 J J ^ t (t ;w) exp [-itanfl p • (u - v)] du dv J (14) 

Since Mj in the large k limit is a function of k only and the total moment M = Mj + M 2 , and since 
by the definition (3) M is symmetric with respect to a k and Ak interchange, then we must have, in 
the limit, 

lim M(k, Ak,«) = Mj(k) + Mj (Ak) (1 5) 

k-*» 

Note that the limit does not have to be taken in order that (IS) hold true, at least approximately. 
The decomposition (15) will be valid in all cases provided that k and Ak are small compared to k 
(condition (3 W )). 


If we defme the surface wave vector 


K s 2xsin0 g 

(16) 

and the specular slope vector 


s = tanfl £ 

(17) 

Then (14) can be written somewhat more compactly as 


Mj(K) = (2r o cos 2 0)” 4 j A(w)exp(-iK* w)dw 


X |(2ir)“ 4 j* j* ^ t (t;w)exp[-B*(u -v)]dudvj 

(18) 


where now the characteristic vector t has the components 

t = (u, -v, Kcotfl, -Kcot0) (19) 

Let us write the solution ( 1 8) as 


Mj (K) = (.2t 0 cos 2 6)~ 4 J A(w)S(K;w)exp(-iK*w)dw (20) 

where 3 represents the factor in braces, viz.. 
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S(K;w)»(2ir) 


f£ 


Kcotfl, -Kcot0; w) 


• expM|'(u-x)l d^dv (21) 

An alternative form of the solution may be obtained by using eq. (13) for ^ and integrating over y 
and v: 

E(K; w) - JJ p^ s, f 2 . J 4 ; w) exptiKcotftf 2 - f 4 )) d^df 4 (22) 

Wc shall find the form (21) more convenient to deal with. 

Kxcept in one particular (to which shortly), the above derivation is essentially identical to Bar* 
rick's ( 1 968b) second-imuncnt derivation. We have used Barrick's method because of the several ap- 
proaches available it is the most direct and to the point. Yet, the reader may feet that the derivation 
lacks rigor; that certain steps require further justification. Rather than attempt this justification here, 
we offer instead in Appendix B an alternative, somewhat more rigorous derivation which sould, by 
analogy to the present derivation, justify it. 

Hie peculiar feature of the fourth-moment solution alluded to (i.c,, the feature of the fourth- 
moment solution that has no analogue in the second-moment problem) is this; the mathematical 
object claimed to be the solution does not exist Die joint pdf P t is singular at the origin w « 0, and 
the integral M, will - gem rally, it is supposed - be divergent. Consider the second moment of back* 
scattered monochromatic power, i.e., consider 

<lS(k)| 4 > ■ M(k. 0, 0) ~ 2M,(0) (23) 

Setting K ■ 0 in ( 1 8) and (22), and in (22) integrating over J 1 and f \ we have 

M,(0)« (2r o eoa l 0)“ 4 J A(w) P vl J (Vf 4 (s.s: w)dw (24) 

where p v( j vf < is the joint pdf of slopes alone. For example, consider a one-dimensional scattering 
situation, where Jtx) is the height profile and J x (x) the slope. If f(x) is a stationary Gaussian process, 
then the joint pdf has a first-order singularity at the origin, i.e., as w-*0 the joint pdf. 


l' f J t 4 (s, s;w) ~ (2 a)"' 4 p f (s) • (g|w|) 

*X*’X *» 


1-1 


(25) 


where g ■ (J 2 X > is the rms curvature. Thus for Gaussian surfaces, the moment M ( diverges logarith- 
mically at the origin and consequently the variance of scattered power is infinite. 

This singular behavior of the solution is a consequence of the point source approximation that 
is implicit in the PO integral formulation (4). The failure of the point source approximation can be 
understood in terms of the statistics of surface cuivature. In the point source approximation (r 0 
the power backscatteied at a specular point is directly proportional to the product of the principal 
radii of curvature, or the same, the reciprocal of Gaussian curvature. Longuet-Higgins (1956) has 
derived the pdf of reciprocal Gaussian curvature g' 1 for normally distributed surfaces and shown 
that the second moment of the distribution is infinite. We surmise that the unboundedness of the 
second moment of g -1 is not a property peculiar to Gaussian surfaces. Rather, we suppose that the 
variance of g” 1 will be infinite for a large class of random surfaces. It can be shown that the singu- 
larity can be removed by including the finite dimensions of the source. When the finite size of the 
source is accounted for, the "singularity" at the origin is so reduced that its contribution to M, (£) 
actually becomes quite small. 

For reference, before concluding here, we write down the result for the second moment 
<|S(k)|* >. From (4), using the same method we have used for the fourth-moment, we have for the 
second-moment in the high-frequency limit 

, /Vcc 4 #\ 

lim <|S(k)|*)» l lp vf (s) (26) 

k-n» \ 4rJ / 

where the illuminated area A 0 ■ /G 2 dx. In terms of the conventional backscatter cross-section, 

o£ d £ f (4irr 2 /A 0 X|S(k)| 2 > (27) 

we have 

o“ ■ trscc 4 0 p vt (s) (28) 

In practice, it is found that the GO limiting form (28) provides a good fit to near-vertical microwave 
backscatter observations provided that some account is taken of the effects of diffraction for finite 


k (Wentz, 1979). Diffraction can be accounted for by treating p vf as though it were die pdf of 
slopes in a smoothed sea surface, in which waves smaller than the diffraction limit (e.g., capillary 
waves) have been removed by a spatial filter (cf. Brown, 1978). Thus, in (28), it is just as well to 
keep the finite frequency tag on o°, as in practice, a k-dependent pdf is involved. 

Following we will perform some calculations in which we assume a Gaussian pdf of orthogonal 
slopes, 

P vf “ (2*)" 1 (detm)-* 4 exp (- ~ m;J s a s fl j (29) 

where mis the covariance matrix of orthogonal slope components, 

a "' 2 

and m" 1 are the elements of the inverse matrix. In (29) the Einstein summation convention is em- 
ployed. A crude fitting of (28) with the Gaussian pdf to the 2. 1 5 cm backscatter data of Jones eteL 
(1977) gives a 2.1 5 cm radar-effective total mean-square surface slope m* that is approximately a 
linear function of surface wind speed U, 

m* * 0.0025 U [ms" 1 ] + 0.01 (30) 

(Ku-band) 

where the slope variance m^ = m aa . The eyeball-fitted result (30) is essentially the same as Wentz’s 
(1979) more thorough analysis of the Jones etal. data. Eq. (30) gives radar-effective slope variances 
in the wind speed range 1 0 to 25 ms" 1 that are approximately 60 percent of the optical values re- 
ported by Cox and Munk ( 1 954). The Ku-band radar-effective slope-variance thus lies between the 
“clean” and “slick” surface observatiowof Cox and Munk (1954). Noteworthy is the agreement of 
the 60 percent figure, here derived from active radar cross-section data, with the percentage slope 
variance inferred from passive radiometric measurements using a GO model of emissivity (Wilheit, 



? ■-» • ' - ' - '■' y i fi " 7^ ■■ n m .in j » " , " i 1 n a >■ /i iji 1 1 
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4. Asymptotic Evaluation of Mi (K) 

We are considering wind-d riven seas in which there is a continuous spectrum of roughness (i.e. 
wave height) from the scale of the dominant wave (MOO m) down to the scale of the diffraction 
limit (~0.1 m). If the incidence angle is not too large, there will always be a large population of 
scatterers (i.e. specularly reflecting wave facets) existing on scales smaller than the dominant ocean 
wavelength. The modulation of the backscatter by the large ocean wave (or spectral ensemble of 
large waves) can then be considered to be in the form of a small perturbation, and an asymptotic 
expansion of £(K; w) about w a *> is appropriate. Let the large-wave steepness, or significant wave 
slope, be defined by 6 0 « KqO where Kq is the wave number of the dominant wave and o is the rms 
wave height. Then the expansion should be accurate if 6 0 is small compared to the total rms slope, 
m Q . This will usually be the case. Thus if fi 0 B 0.0S (a value characteristic offully developed seas) 
and m Q * 0.19 (wind speed of 10 ms -1 by (30)), then 6 0 /m 0 » 0.26. 

M, (K) (eq. 20) is expanded by first expanding E(K; w) (os given by eq, 21 ) and then term-by- 
term Fourier transforming with respect to w. £ is expanded according to Longuet-Higgins’ (1963) 
method for generating the non-Gaussian pdfs of variables in a non-linear Stokesian sea. Indeed, we 
have chosen the form (21 ) for £ in order that we could parallel Longuet-Higgins’ (1963) development 
And so lay the foundation for a theoretical model of Mj (K) capable of entertaining non-Gaussian 
wave statistics in a rational manner. 

As can be seen from (13) the coefficients in a Taylor series expansion of the characteristic func- 
tion correspond to the moments of the distribution, i.e., if 

i i 2 

(1; w) « 1 + — n\ l \ + — gjj 1 * tjtj + . . . (31) 

then the jj coefficient of order n is the n-th order moment, 

\ { - )m f ( V* V p i ( * ; - )d * 


(32) 
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The cumulants of the distribution, by definition, correspond to the c coefficients in the expansion of 

C t (|;w)-tn# t (i;w) 

"TT c * l>t| ' 1 ’ *1? c P t ‘ t i + -*- (33) 

Longuet-Higgins (1963) shows how, in general, it is easier to calculate the cumulants (as opposed to 
the moments) from the non-linear, Eulerian water wave equations. Further, an advantage in working 
with the cumulants is that if £ is nearly normally distributed, then the third and higher-order teams 
in the C series will be small. (This is so since for | exactly normal, C is given exactly by the first two 
terms of the expansion (33).) Considerable simplification of the problem is possible if we restrict our 
attention to statistics of the third-order or less. Thus, since J (cf. eq. 1 2a) has aero mean, c (1 > « #|W ■ 0. 
it follows that the second and third order moments and cumulants are identical, as can be verified by 
a term-by-term comparison of (31 ) and the expansion of ^ ■ exp(Cj). Thus, if we restrict ourselves 
to third-order statistics, can be approximated by 

« cxp(Cj) - exp 

X ll * y, (J4) 

All elements of the covariance matrix “ <$, (j > can be expressed in terms of the wave height co- 

variance function R(w), 

R(w) J « >-<£(*•► w>f(*)> (35) 



or various derivatives thereof, as can be verified, for example, using methods outlined in Papoulis 


(1965, p 314 AT): 



( 36 ) 
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where we have defined o i ■ R(§), R ■ d R(w)/9w Q etc. and m afl ■ -R a(J (2), and where a and 0 
range from 1 to 2. 

Let where d and o denote the diagonal and off-diagonal blocks of Since 

■ 0(mj) and ■ 0(6 2 ) for w removed from die origin, i.e. on the order of wavelength, we 
can then expand (34) for small 

~ exp (■ \ x [■ - \ < H 2) Vj 

■ $ ^(jk Wk + 0#i kV • * J ^7) 

where and are functions of w, but is not. Similarly to (36) we expect that the ele- 
ments of are all derivable from a knowledge of the third-order mean lagged products of surface 
height Q(r, s) ■ <f(x + r) f (x + s)f (x)>. Thus, die FT of g* 3 *(w) should be expressible as various 
operations on the bispectrum, the FT of Q(r,s) (Hasselmann etal. (1963). In principle, the bi- 
spectrum (and hence the elements of g* 3 ) and the FT of g r3) ) can be derived from the height of 
spectrum F(K) = FT{R(w)} using a second-order Stokes type expansion of the nonlinear equations 
for free gravity waves according to the theories of Hasselmann etaL (1963) and Longuet-Higgins 
(1963). 

We are thus proposing here that the appropriate model for the near-vertical wave-spectrum 
measurement problem is one that (a) models scattering according to GO and (b) models the sea- 
surface accordi.ig to the Longuet-Higgins ( 1 963) and Hasselmann et al. ( 1 963) theories. That this is 
a realistic model for near vertical backscatter is indicated, for one, by the success of such a model in 
predicting the average impulse response of the sea surface to 3 cm radiation at vertical incidence 
(Jackson, 1979). Moreover, regarding the appropriateness of the GO approximation: if there is 
reason to suppose that GO is a good approximation in the average impulse response problem, then 
there is even more reason to suppose that GO is a good approximation in the present problem, i.e. 
that of modelling the spectrum of the impulse response. This is because in the present problem we 
are only interested in the variation of the backscatter over the large wave profile, and not in the 
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absolute power level. Thus, while the diffraction component may be small compared to the GO com- 
ponent (cf. Brown, 1978), the variation of the diffraction component will be smaller still compared 
to the GO component. This is because the diffracted fields are very diffuse, and can therefore con- 
tribute only negligibly to the modulation via the geometrical tilting mechanism. 

We shall be neglecting the non-Gaussian terms in although there is noteason to suppose 
a priori that the g* 3 * terms arc negligible compared to the second-order terms in Thus, for 
example, consider small values of tan0/m o . Then the terms such as R(w)K 2 cot 2 0 predominate. 
Thus the linear term is O(K 2 o 2 cot 2 0) * Offi^cot 2 #) and the second order °ji^ 2 ) term is O(8£cot 4 0). 
The corresponding gP* term will be O(Xo 3 K 3 cot 3 0) ■ O(X5 3 cot 3 0) where X * <f 3 )/o 3 is the height 
skewness coefficient. In general, X and 6 0 are proportional to each other (Huang, 1980; Jackson, 

1 979). From Longuet-Higgins’ ( 1 963) theory, if one assumes a Phillips spectral form, one gets 
X * 48 0 . This relationship follows from Jackson ( 1 979) if an apparent error of a factor of two in 
the second-order height profile £ 2 is corrected for; i.e. if «- f 2 + 2 (D. E. Barrick, personal com- 
munication, 1980). The ratio of | j* 3 * terms to second-order °g( 2 ) terms is then of the order of 
X/6 o cot0 = 4tan0 * 0.7 1 for 0 “ 10°. Hence, a priori, the neglect of the g* 3 * terms is not justified. 

Proceeding nevertheless with a Gaussian surface model, using (36) and (19) (for the definition 
of t ), we have 

4> t ~ V'* vr (u) (-v) exp(-(Kocot0) 2 1 

X (l +(Kcot0) 2 R + Kcot0R a (u a -v a ) 

-R 0<) u 0 v fl + second-order terms in < V 23 } (38) 

where = expi-iim^u^). The-sccond-order terms are carried in Appendix C. Now making 

use of the identities, 


P f (s) « — j* * (u)e^**B du 

(2v)‘ J 


3p vr > f , . . -*5 -S 

3 *. J U ”’ tvI ' r dl ' 


( 39 ) 
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and so forth, we have for S (eq. 21 ) 


S ~ p 2 f ( 5 ) exp (-(Kocotfl) 2 1 • |l + 

, p .« 

+ (Kcotfl ) 2 R + 2iKcot0 -j- R # 

P.P, A n\ I 

R + second-order terms in °u w 1 

p2 * I 


where for short in the series part of (40) we have let p * p vf (s) and p a a 3p/9s Q , etc. The second- 
order terms are carried in Appendix C 

The moment (K) as given by (20) is the FT of S as seen through the lag window A(w). If 
A(w) is broad compared to the correlation scale of the waves (as will always be the case in satellite 
measurements), then the effect of the finite window on the modulation spectrum can be igonred and 
we can write as 

Mj (K) = (jr/2r 2 cos 4 0 ) 2 p 2 f (s)exp {-(Kocotfl) 2 1 

x [FT {A(w)} + A(0)* FT {linear + higher-order terms in (40)}] (41) 


where the two-dimensional Fourier transform, 


FT 


def. 1 


i j*{. . e -i £* 2 dw 


Let the reference axes x j and x 2 be fixed in the antenna beam; that is let x j be the surface range 
coordinate, parallel to g and K, and let x 2 be the orthogonal azimuthal coordinate. Assume a Gaus- 


sian gain pattern. 


G(x) • exp [- \ (Xj/Lp ) 2 - -y < W*] 


then (using ( 1 1 ) for A(w)), 


FT {A(w)} *(2a )" 2 A(Wj , w 2 )exp(-iKwj )dwj dw 2 

»(A 0 /2)exp[-i(KL p ) 2 ] 
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Then, malting use of the second-moment result (25), noting that A(Q) ■ Aq/ 2 ■ »L p L p /2, we write 
(41) in the form 


M 


,(K)» ~ I<|S(k)l 2 >l 2 j^=e-“« Icl *<s> a ♦P luod (K) 


(45) 


where the modulation spectrum, or spectrum of the impulse response, is 

P mod (K) ■ (v/T»F/L p )expl-(Kocot0) 2 ] • FT (linear + higher-order terms in (40)} (46) 

The directional height spectrum is defined by 


F(K) d - ’FT (R(w)} 


(47) 


Then, since FT (R J « iK a F(K) and FT{R >ofl }«-K a K fl F(K), the expansion for P rood becomes 


mod 


(K) = 


y/lv 


.-(Kocotfl) 


'It 


P« K a 
cot 2 0 - 2cot0 — • — 
P K 


P 2 K 2 


K 2 F(K)+... 


or 


P mod (K)» y* e-< KoCoW)J • |(cot0- ^?)*K 2 F(K) + ... 


(48) 


(49) 


where dKnp/ds = (K a / K) * (p a /p) = g • V,p. The second-order term in the expansion of P mod is 
given in Appendix C, 

Some features of the solution are noted. If the exponential factor (the square of ) is 
neglected (which neglect is consistent with the neglect of the second-order terms) then the linear 
ized solution (49) yields the same result as the simple tilting model (Fig. 2). The dimensional pro- 
portionality factor L~ l relating P mod (K) * m to K 2 F(K)& m 2 is the reciprocal of the azimuthal 
beam spot dimension. An L” 1 dependence will exist whenever L p exceeds the lateral correlation 
length scale (e g. crest length) of the surface waves. Large L p means small signal. This loss of modu- 
lation signal power with increasing L p can bt thought of as the price paid ;ar high directional reso- 
lution: of the ensemble of Fourier component surface waves contributing to a total modulation 
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power, the radar is isolating only a small subset that are travellihg in a small direction band. The 
dependence can also be interpreted as follows: in a real short-crested sea, the waves are running out- 
of-step. Hence, modulation signal power adds roncoherently across the beamwidth resulting in an 
dependence. The total signal power on die other hand adds directly as L # . Hence the relative 
modulation power goes as LJ* 4 , and die spectrum of power as L" 1 . 

To gain some feel for the effect of the second-order terms, sample calculations are carried out 
for a simplified two-dimensional scattering situation using eqs. 49 and CS of Appendix C. The results 
shown in Fig. S indicate that there is a fairly good measurement fidelity to the wave slope spectrum 
in the incidence angle range 8° to 1 5° for wind speeds > S ms** 1 and significant wave slopes 8 0 < 0.1. 


5. SNR Analysis: The Modulation Noncoherent Case 
Narrowband Dual-Frequency (DF) Technique 

The appropriate signal-to-noise ratio (SNR) in the narrowband DF technique is the DF correla- 
tion coefficient (cf. Barrick, 1972; Plant, 1977). Thus, from (3), (IS), and (4S), we have 


SNR def._M(lc 1 K^ _ j 
I<|S(k)| 2 >1 2 

= (v^T/L p )P mod (K) ('**; 

In (SO), note that P mod is the one-sided spectrum, i.e., we have let P mod «- 2P mod . Here, and in die 
SNR calculations that follow, we use only the linear part of P mod with the exponential set equal to 
unity. (This assumes that 5 o cot0 « 1 , and is consistent with the neglect of the second-order terms.) 
Let the one-sided spectrum F(K) be given by the Phillips spectral form with a cos 4 4> spreading 
factor. 


F(K) ■ (0.005) (8/3jt)cos 4 <1 > K“ 4 , 


K>K 0 ~g/U 2 

l*l<»/2 


(51) 


[ 0 otherwise 

This spectrum defines a sea with £ steepness 6 0 * K 0 o s K 0 (/F(K)dK] ** * 0.0S. For upwind (up- 
wave) looks 4* ■ 0, and F(K) * 0 0042K" 4 . From (28) with s ■ tan0( 1 ,0), we have dCnp/ds * 
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-tan0/m| O ; hence, 

P mod (K > * (v^7/L # ) (cot® + tan0/m, o ) 2 (0.0042)K’ 2 (52) 

Let K * Kq ■ 2w/200 ra, and corresponding to which, U * Vg/Ko a 18 mi* 1 . If we use the Cox and 
Munk (1955) upwind/crosswind slope variance ratio m Q | a* .66 m t0 , then m (0 * m|/l.66. From 
(30), we have m| at 0.055 and so m, 0 at 0.033. Let 6 ■ 1 2.5° and L^ ■ 2.5 km; then 

P mod 31 («0 -3 nT 1 ) (4.5 + 6.7) 2 (43 m 2 ) - 0.54 m (53) 

If ■ 5 km, then from (50), the SNR B -36 dB. In addition to the SNR, defined here as the ratio 
of the modulation signal spectrum to the square of the average power (equals the self-clutter, or ran- 
dom power fluctuation ariance), we need a figure to express the signal strength relative tp the aver- 
age backscattered power. Call this figure the contrast ratio, 

CRd “^ P mod^ AK ) ,4 "^ NR ) V4 " , - 6% (54) 

where the last follows since AK = >/Jirl L p . 

Short Pulse (SP) Technique 

The frequency-domain analysis of the SP system of Fig. 3a is carried out in Appendix A. In the 
transformation between frequency and wave number domains, we have the identities p ■ kc, o> ■ kc, 
and Av * A kc. Since K « 2«csin0, etc., the moment result, eqs. (15) a::d (45) becomes in the fre- 
quency domain 

M(v,Av,<a>) ■ M,(^,<*>) + M,(v,Av) (55) 

where 

M, «(^/T r )I<lS(v)| 2 >l 2 {6(w) + P m ^(w,*)} (56) 

and wiiere T r = (2L p /c)sin0 and where for large T r we have let 

5(w) ~ (T r //2ir) exp — (wT r ) 2 J 
where 6(u>) is the Dirac delta function. 


24 


The backscattcr of a short pulse (i.e. “short" compared to the ocean wavelength) consists of a 
weak modulation m(r) riding in a broadband, exponential “selfclutter" process w c (r) of bandwidth 
equal to the pulse bandwidth dp- (The exponential statistics of w c (r) follow from an assumption of 
normal statistics for the unmodulated complex field.) Thus, the backscattcr power is of the form 
p(r) ■pll+ m(r)) w c (r). Using either this simple time-domain model, or using die moment analysis 
results of Appendix A, one finds that the spectrum of p(r) has the form (e,g., assuming a Gaussian- 
shaped pulse): 

Pfw)« ^KlSMW e-*<«/'»> 2 |6(w) + ~rj- ♦P mod (u.,*)| (57) 

where the approximation is valid for small modulation, i.e. for / P mod ( Av) Av « 1 . Eq. (57) is 
sketched in Figure 6. 

Note that, if we let d p -* then the transmitted pulse becomes an impulse. Thus P mod is just 
the normalized spectrum of the surface impulse response. The SNR appropriate to the SP technique 
is the ratio of the signal spectrum to the self-clutter noise spectrum, i.e. the SNR ■ dpP mod . 
Now, some or all of dp may be used for range resolution. A minimum pulse bandwidth is required 
to achieve a desired range resolution (e.g. 25 m on the surface). Bandwidth in excess of this mini- 
mum serves to reduce the level of the self-clutter spectrum. If we define an equivalent range resolu- 
tion cell A, such that if all the bandwidth were utilized for range resolution, i.e. A p s c/2d p sin0, then 
we can express the SNR in a form similar to (50). If we take as ac example Ap • 5 m, then (again 
using (53) for the estimate of r mo< j) 

SNR - (%/27/Ap) P rood (K) - -6 dB (58) 

which is 30 dB above die narrowband DF SNR, i.e. grater by the factor L p /Ap. The contrast ratio 
is defined as the rado of the rms signal level to the dc level; from ( 52) we have 


CR« 



1 mod'— 


(K)dK 


It 


13% 


(59) 


v 
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As£hc SNA can always be impro'^by poise integration (as we shall shortly show), the CR could 
turn out to be the more important figure in design considerations. The 1 3 percent figure is quite 
good; what is more, the CR will still be quite good for shorter waves since CR <* Thus, e.g., 
for a fully developed (5 0 * 0.05) sea wave of 50 m in length, the CR is reduced by only a factor of 
two from the 1 3 percent CR of the 200 m water wave we have used as an example. 

Wideband DF Technique 

The wideband OF technique can employ a variety of waveforms including SP waveforms. 

(Thu*, if we transmit pulses shorter than the dominant waves then, necessarily, all the required dif- 
ference frequencies are carried by the pulse.) If the transmitted waveform FT is, for example, assumed 
to be constant over the bandpasses of the filters Hj and Hj of Figure 3b then, as shown in Appendix A, 
the strength of the DF signal at the output of the correlator will be governed by the combined band- 
width of the Hj and Hj analyzing filters in each channel (cf. Figure 3b). Owing to the action of 
the dc blocking filters and the fact that the self-clutter fluctuation power is uncompleted between 
frequency bands Av > 2ff/T r , the expected value of the output of this DF system will contain neither 
dc nor self-clutter components. 

So. what is the appropriate SNR for this system? And how should tills system's performance 
be compared to the SP system? Clearly, we have to look at the detection problem more carefully. 
Parpen and Shiren ( 1 956) went some distance toward this by establishing a detection criterion for 
the SP system. Unfortunately, they stopped short of a like analysis of ihe DF system. Lacking a 
rigorous analysis by Parzen and Shiren, we will do with some back-of-the-envelopc calculations. 

Assume that the signal levels arc within Uie dynamic range and resolution capabilities of tnc instru- 
ments. .Assume further that we can ignore the sampling variability in the estimate of P mod compared 
to the variability in the self-clutter spectrum. Then, successful measurement depends only on the 
strength of the signal relative to the sclfeluttcr fluctuation levels about thetr expected values. Con- 
sider the SP system. The signal at a frequency u/j appearing through a finite analysis window of 
bandwidth is ^ 1| P tT1VH) (co ) ). Tlie mean noise level is /3 41 /V3ir0p. The noise standard deviation in 
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| a measurement of total time-bandwidth product TBP = N p T r /3 a{ /2ir where N p is die total number of 

independent pulses, is (y/7wfi p )~ l (TBP)-* 4 . Hence, if we define the detectability D* as die ratio of 
the signal to the noise standard deviation, then 

i 

In the DF system, the signal out is proportional to 0 ai jP mod (Ai>y). Assume that the mean square 
fluctuation of cross-correlated power as a percentage of the product of mean channel powers is simply 
TBP-* 4 where TBP = N p /J a yT r /2ir. If this assumption is correct then the D* for the DF system is 

Dfo, = (N p T r )»S$ P mod (Av*> ««) 

If die DF analysis bandwidth 0 a y is opened wide (i.e., made comparable to the Avy), then D^p is 
comparable to the Dg P . Any further analysis of SP and DF techniques must consider thermal noise, 
antenna gain and other tradeoffs. This is beyond the scope of the present work and is, after all, 
academic, considering that the SNR can be improved immediately by making the detection coherent 
with respect to the modulation signal. 

6. A-H Filtering: The Modulation Coherent Case 
Alpers and Hasselmann (1978) showed how by filtering the complex DF signal 

D(fc,t) » S(k,t)S*(k - x.t) (62) 

at the appropriate Doppler frequency fl = K • V, a dramatic increase in the narrowband DF SNR 
could be had. We show how, if the same type of processing is applied to an SP system, the same 
dramatic increase in SNR results. As the narrowband DF results follow immediately from the more 
general moment results in the analysis of the wideband systems, we proceed directly with the wide- 
band analysis problem. 

» 

I Consider the SP system of Figure 4. The output of the square-law detector is p(r,t) = |e g (r,t)| 2 . 

i 

The FT of p with respect to fast time is denoted P(w,t), 

E g (i»,t)E* (v - w,t)dv 

! 


Aw,t) = 


-1 


(63) 
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where by (2) E g (»\ t) * S(M)E o («0. The pulse spectrum is defined as P(w) » <|P(w,t)| 2 >. The simple 
filter in Figure 4 functions to integrate pulse returns in surface-fixed range bins. The output of tire 
filter after an integration time T int is 

I f Tjnt , 

P(f)“qp— I p(r + Tt,t)dt (64) 

*lnt Jo 

where r * 2r /c. The spectrum of the filtered pulse return is P(w)“ <|?(«,t)|*> where P(ox)> is the 
FT of p (r). Hie spectrum of the filtered pulse return can be expressed as 




w(n-K*v)P(w,n)da 


where the filter window W(ft) 


W(ft) = 


sin(ftT int /2) 

(»T lnt /2) 


»5) 


( 66 ) 


and P(u>,ft) is the post-detection pulse- Doppler spectrum defined by 


P(w,ft) = 



<P(u>,t) P*( w,t + At)>e in A 1 dAt 


(67) 


The Dopplor frequency ft ® K • V = -our. If 4> = 0 coincides with the velocity vector as in Figure 1 , 
then 

ft = K-V = KVcos4> (68) 

If the generalized moment is defined as 

M (v, t>\ w; A t) = <S(v, t) S* (i> - w, t) S* (»>', t + A t) S(v' - w, t + A t)> (69) 

and the FT of M with respect to At, 


1 

N(v,v’,w;ft)® — I M(*>,»»',w: At) e lnM dAt 
2 IT J_oo 


( 70 ) 
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then 


m 

p( Wl n)- JJ N(*y , u; n> ^ w e* - w) e* E 0 y - w) d»> d*' an 


Now let us again equate the fast frequencies *\ v\ and u with die wavenumbers k, k', and k, (by 
¥ ■ kc, etc.), so that in terms of wavenumber M ■ M(k, k', k) ® M(k, Ak, k). We calculate M with 
S(k.t) given by PO integral solution (4) in the Fresnel approximation to the phase (5). where we 
neglect the term in height. kj 2 /r„. Proceeding exactly as we did before in the fareone approxi- 
mation, we form the moment as the fourfold integral over the dummy space variables & 1 . . . and 
interchange expectation and integration operations. Again we recognize that for large k, significant 
contributions to the moment are made only in the vicinity of the “stationary** points defined by (8). 
Hie moment M is thus decomposable as M * Mj + M 2 where M t and M 2 represent the contributions 
from S, and 5,. Again, we consider die integration over die volume surrounding the S, points. 
Transforming to thr lag variables u, v, and w as per (9), and taking the limit, we get for the contribu- 
tion from S, 

lim M, ■ M,(K; At) ■ (^cos 2 #) -4 j f G 2 (x + w)G 2 (x - V At) 

* * 

• S(K; w)expHK*(w + VAt)l exp{i(x/r 0 ) 

• (cos 2 0(x, + Wj) 2 -(x, -V,At) 2 +(x 2 + w 2 ) 2 

-(x, -V 2 At) 2 l}dxdw (72) 

where we have taken the coordinates x, and x 2 of x to be radar-fixed coordinates (x { in the plane of 
incidence). Let G be the separable Gaussian pattern given by (43 V The integrating over x we get 

M|(K;At)«(2r e cos 2 *r 4 J B(K;w+VAt) 

• 5(K; w)cxp|-iK*(w + VAt)] dw (73) 

where 

B(K; w) * (Aq/ 2) exp ~ b 2 w 2 j (74) 


» 


and 

bf * b 2 « L“ 2 + (KL p /2r o sin0) 2 co$ 4 0 

0 sin») 2 

^■ir^L^/2 (74*) 

and w icre again the summation convention is applied to repeated indices. 

Since M ■ M, + M 2> then N as given by (70) is ■'imilarly composed, i.e., N * Nj + N 2 where Nj 
and N 2 are the respective PTs of M t and M 2 . Ftom (70) we have 

Nj(K; fl) « (2r 0 cos 2 ®)' 4 Jdw S(K; w) exp(-iK* w) 




(2ff)-> dAtB(K;w+VAt)exp(i(n - K*V)Atl 


(75) 


Let W » (bj w, , b 2 w 2 ), V « (b, V, , b 2 V 2 ) and - K* V. Then integrating over At we get 

(Ao/2) f 

Nj(K;Si)» — I dw3(K; w)exp(-iK*w) 

(2r o cos 2 0) 4 J 


exp 


[-4^./^] r , , 

Tliy— ■' xp l‘2* staJ(t - ) J 


• expUW'cos(K. W)ajV] 

From (70) and the inverse FT relationship, it follows that 


(76) 


M,(K;0)- 




J* N,(K; I 

Crocos 1 ^)- 4 j B(K;w)2tK;w)c-*K!? dw 


(77) 


which is the same as the Ftaunhofer solution (20), except that the Fresnel lag window B(K; w) re- 
places the Fraunhofer lag window A(w). Since B is not significantly sharper than A and the window 
can still be considered to be broad compared to the scale of variation in 2 (i.e., broad compared to 
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the dominant wavelength 2w/K 0 ), then M| (K; 0) can be regarded as identical to the Fraunhofer 
moment M|(K). (Note that B(0; w) ■ A(w) and B(K;g)» A(Q) * Aq/ 2.) If the detailed behavior 
of N^K; ft) near the origin K a 0is ignored, then it is apparent from (76) and (77) that Nj is well 
approximated by 


N^K;Q)«M 1 (S> 


exp [- i (£l,/^) 2 ] 


where M| is given by (46) 4 tseq., and 

fl m - f * iKl»(b ft V a )* “(bjcos^ + bjsin 2 *)*^ ° 9) 

where 4> is the radar azimuth relative to the velocity vector V. In A-H, only the lino-broadening due 
to wave-front curvature was considered; i.e., it was assumed that fi m =* «LV/r 0 . However, to be ac- 
curate, the broadening due to the finite footprint must be accounted for also. Thus, compare V/L 
and KLV/r 0 for typical aircraft and satellite geometries and ocean wavelengths: both terms are 
comparable. 

The moment M(k, Ak, k; At) is not symmetric with respect to a k and Ak interchange, and so 
the M 2 component (deriving from the S 2 stationary points) cannot simply be equated to Mj, as was 
possible in the time-independent case. Rather than labor through a solution for M 2 in these pages, 
we merely outline the key steps in the solution and indicate the appropriate assumptions and approx- 
imations: lag variables appropriate to the S 2 integration are defined, and the k-+«° limit is taken in 
the formation of S; on account of the moderate bandwidth assumption (3"), the Fresnel phase terms 
in k and Ak are neglected; assuming a large footprint area, we neglect the variation of S near the ori- 
gin and let E(0,w) at £(0,<*) a p 2 f (s) in the integration over w. Then, on taking the FT of M 2 , we 
get 


N 2 * l<|S(k )| 2 » 2 exp [- ^ (qAKL p ) 2 ] 

exp (syu d ) 2 ] 


sfT* ft, 


txpdcn j&l) 


(80) 
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where 


q * |sin$|‘ t(L p /L^) 2 cos 2 $ ♦sin^]"* 4 
0 d ■ (2kV/r 0 )* [L 2 cos 4 0cos 2 $ + L 2 sin 2 4>] % 
e a (2kV/r 0 )* AKL 2 cos 2 0cos4 (81 

and where AK a 2Atam0 3 2(Av/c)sin0, and where again H 0 = $2 - K* V. Asa check on (80) we 


J N 2 dSl»l<|S(k)| 2 » 2 expJ-l(AKL p ) 2 ] 


which is Mj(AK) with the weak modulation contribution neglected. 


7. Results: Narrowband DF Technique 


The spectrum of the narrowband DF signal, eq. (62), is simply N with AK set equal to zero 1 in 
Nj. If a signal of finite record length T int is passed to a spectrum analyzer or, equivalently, if a 

A 

Bartlett lag window is applied to a measured covariance function M in a Tukey analysis (Blackmann 
and Tukey, 1 9S8), then the measured spectrum will be (on the average) 


P(K) = Jw(n-K*V)N(k,K,0;Sl)dSl (8 

where the spectral window is given by (66). If T int is chosen such that <Wr in t « 2jt, then from 
(78), (45), and (80) we have approximately 

P(K) s p, + P 2 « W(0) Jn, dfl + N 2 (0) J Wdft 

For example, let us set T int = /3L* . Then the measurement signal-to-noise ratio is 


SNR.?,/?, ^ P mod (K) 


where we have let P mod become the one-sided spectrum; i.e., we have let P mod *- 2P mod . With the 
same parameter values used previously in the modulation noncoherent analysis, i.e., 6 = 1 2.5°, 
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L p « 5.0 km, ■ 2.5 km, |K| ” 2*/2QO m, U ■ 18 ms” 1 , we have as before (>/I»/L p )P mod (K) ■ 
-36 dB. Now, if k “ 2w/2 cm, r 0 ■ 700 km and V s 7 kms -1 , then e.g. for broadside looks, d>“±90°, 
we have 

P m a KV/L p ) 2 + (KL,V/2r o sin0) 2 l* 

» 17.84 s" 2 + 3.29 s“ 2 J % ■ 3.34 s” 1 
U d » 2kL 0 V/r o - 1 .57 X 1 0 4 s (86) 

The filtering gain is (2tr r 3 ^(0 d /0i«) * 25 dB. Hence, the final measurement SNR is 

SNR DF"‘ 36dB + 25dB»-U dB (87) 

Of course, this SNR figure can be improved upon by designing a spectral analysis having a better 
“match" to the resonance peak (e.g., by using a matched filter). However, it is apparent that, at 
best, given reasonable parameter values, the SNR will not much exceed 0 dB. 


8. Results: SP Technique 

Making use of the transformations (55) et seq., we express our results (78) and (80) in terms of 
frequency. For |<o I > 0 (i.e., ignoring the dc), we have N(v, Av, to; N, + N 2 , 


Nj »KlS(v)| 2 >l 2 exp|--i^^ J • “^P mo <i(w.*) 

N 2 a l<|S(v )| 2 >] 2 expf- 1 CX P l-y (qAvT r ) 2 J t\piieSlJP A ) 


(88a) 

(88b) 


If the pulse bandwidth is large compared to (qT r )'* , then N 2 behaves as a delta function in Av= v-v'. 
Hence, since <|S(v)| 2 > is practically a constant over the ;»u^e bandwidth, eq. (71) becomes 

P(w, O) = P, + P 2 

■N,-|/E 0 (i»)Ej(i>-w)di‘t* 

+ 1 / E 2 tv)E^(»»- w)dvl •/ N 2 dAv 


( 89 ) 
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Making use of the definitions (81 ) and the identity AKL p b AvT, we have for the N 2 integration, 

C exp [“ T 

J N 2 dAv*(v^?/T r )l<|S(v)| 2 >l 2 “ C 


where (from 81), 


q^d “ (2kV/r 0 ) L IsinC* | 


which is »he Doppler spread due to azimuthal aspect variation alone. 

Pulse waveforms are often Gaussian in shape; let E 0 (i>) at baseband be given by 


E 0 (v) ' 


[4 ( -v] 


Vi* 


Then, computing the integrals in (87), we have 


P(w,n) = P 1 +P 2 = — 1<|SC*^>| 2 >1 2 • 

*r 






cxpf-^tn^) 2 ] 






exp - -j (Syq0 d ) 2 ] 
2tr q/3 p <3 d 


The post -detection “pulse-Doppler" spectrum P(w. SI) is sketched in Figure 7. 

Except for looks directly forward (4> = 0°) or aft (<t> = 1 80°). q0j » Hence, in the 
vicinity of the resonance peak at £2 = K • V we can set = 0 in P 2 . Again, if we assume short in- 
tegration times, T int , such that the resonance peak is entirely covered by the spectral window W, 
i.e. ii fJ_ g Tj nt « 2ir, the filtered pulse ••pec trum becomes (cf. 65), 


P(w) = P, + P 2 =* W(0) PjdSi +P 2 I WdO 


Again, assuming T in , = 0]^,. from (93) we have 


SNR sp = P,/P 2 -«(2» r)“ 


J P r mod l 


Again, we assume the same familiar parameter values (cf. Table I). Then, as before, we have 

0 p P mod (w. *) ■ (>/37 /A p ) P mod (K) - -6 dB (96) 

For broadside looks (4» = ± 90°), we have ■ 3.34 s" 1 (from 86a) and q/3 d » 1.57 X 10* s” 1 (from 
86 and 91). Hence, 

SNR sp » 25 dB - 6 dB a +19 dB (97a) 

(broadside) 

For looks ± 1° of forward or aft, 0 m = 2.23 s _1 and q0 d * 2.74 X 10 2 s" 1 ; in this case, 

SNR S p * 9dB-6dB = +3dB (97b) 

(fore/aft) 

9. Discussion 

In the foregoing analysis, the backscattered signals were assumed to be continuous functions of 
slow time t. Such an assumption will be valid only if the pulse repetition frequency (PRF) is greater 
than or equal to the rate of the fastest signal fluctuation determined by the Doppler spread (i.e., de- 
termined by 0 d in the narrowband DF measurement (beam-limited illumination), or determined by 
q/3 d in the SP measurement (pulse-limited illumination)). If the PRF is lower than the Doppler fad- 
ing rate, then signal fluctuation power at frequencies greater than one-half the PRF (the Nyquist fre- 
quency) must appear in the frequency band 0 - Vi PRF. The pulse Doppler spectrum P(w, £2) at 
multiples of the Nyquist frequency, will be folded over onto the interval 0 - V4 PRI (.Blackman and 
Tukey, 1958). This folding over will result in an increased measured signal fading spectrum, the in- 
crease being roughly in proportion to the ratio of the Doppler spread to the PRF. Clearly, in the 
case of the narrowband DF measurement, where the SNR is already poor in the continuously 
sampled case, one cannot very well afford to undersample. In the SP measurement, on the other 
hand, one starts with a generous SNR in the continuously sampled case, and one can easily afford to 
lower the PRF. Clearly, the designer of an SP instrument will be able to select and trade off pulse 
bandwidth, integration time, and PRF parameters with a fair degree of freedom; he will also have 
some flexibility in terms of antenna gain and azimuth scan rate parameters. 
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We have not undertaken to describe or analyze a possible modulation coherent wideband DF 
system. Our experience with the modulation noncoherent wideband DF system, however, would 
indicate that such a system would perform comparably to an SP system. 

We have not considered thermal noise in our analysis, mainly for the reason that thermal noise 
will not be a limiting factor unless it equals or exceeds the average backscattered power level. Avail- 
able transmitters (e.g., the one in the Seasat-1 altimeter (Townsend, 1980)) have enough power to 
deliver a unity-or-better signal-to-thermal noise ratio when modest gain antennas are used and the 
nadir angle is about 10°. 

We have not given any Figures for directional resolution; some actual numbefs may be enlight- 
ening. These numbers can be derived from the Fresnel lag window B(K; w), eq. (74). If we consider 
only the linear term in S(K; w) (in which the exponential factor is neglected), and if we assume 
K » h a , a - 1,2, then it is apparent that the spectral window through which the slope spectrum 
K 2 F(K) is seen is given approximately by the FT of B(K; w) with respect to w. The one-sigma val- 
ues of the Gaussian-shaped window in range (a = 1 ) and azimuth (a = 2) are respectively given by b a . 
Using the familiar parameter values of Table I, and defining directional resolution as the half-power 
window width in azimuth, we have 

A4»=2>^2ftn2b 0 /K 

* 2V2taT |(KL # r 2 + (L*/2r o sin0) 2 l V4 
= (2.36) • (1.6 X 10- 4 + 1.3X KT 5 ]* 

= 0.031 rad = 1.8 deg. (98) 

For the main reason that the diffraction fields in near-vertical backscatter are both small and 
diffuse compared to the GO (geometrical optics) fields. 1 believe that the GO solution for the modu- 
lation spectrum given here (or the same, the non-dc portion of the spectrum of the impulse response) 
is an accurate one, at least for incidence angles less than 1 5°. The expansion of Z(K; w) to the first 
order in non-Gaussian statistics (i.e , 'bird order moments) and to the second order in Gaussian 
statistics is, 1 believe, a reasonable and consistent approach to take. By expanding E (asymptotically. 
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for large lags) one has immediately a linear term thht is directly proportional to the directional wave 
slope spectrum, with a proportionality factor that agrees with the simplest tilting model of reflectivity 
modulation (cf. Fig. 2). Also, by using an expansion technique, one has a method, in principle at 
least, for solving the inverse problem : the raw observations can be used to give a first estimate of 
the true directional height or slope spectrum which, in turn, can be used in a forward prediction of 
the higher-order terms. Since those terms are small, one iteration would probably suffice. Still, a 
means of handling the non-Gaussian terms in a theoretical sense and in a practical sense (i.e., in terms 
of developing a viable algorithm) needs to be invented. 

While the expansion of S(K; w) is certainly valid for large lags w, it does not follow that the 
expansion of M t (K) is similarly valid, as stated in the abstract of this paper. Since S is not accurately 
represented for small lags, then M,(K), as the FT of S with respect to w, will be uncertain in terms 
of the absolute degree of whitening caused by the fine structure of 3 near the origin. Until this 
representation problem has been solved, or at least until quantitative error bounds have been estab- 
lished, the sample calculations given in Fig. 5 must be regarded as semi-quantitative. 

10. Conclusion 

We have shown, theoretically, how rather simple, scanning-beam, microwave radars can be used 
to measure directional dcean wave spectra from satellite platforms. The theory indicates that: 

(1 ) The short-pulse (SP) approach to the measurement is superior to the narrowband dual-frequency 
(DF) approach in terms of measurement signal-to-noise ratio (SNR) and contrast ratio (CR), and in 
terms of the added flexibility in overall design that the greater SNR and CR figures of the SP approach 
afford. (2) Conceiv bly, wideband DF systems may be designed that perform comparably to wide- 
band SP systems. (3) The radar-observed surface spectrum bears rather good fidelity to the direc- 
tional slope spectrum in the incidence angle range 8° < 6 < 1 5 J provided that the wind speed is not 
to low (U > 5 ms -1 ) and provided that the large wave steepness is not too great (6 0 < 0.1). In 
principle, based on the theoretical model presented here, algorithms can be written to remove the 
small amount of harmonic distortion inherent in the specular scatter measurement of the sea slope 
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spectrum. Such algorithms, it is believed, will have to take account of the effects of non-Gaussian 
water wave statistics, as it appears that the effects of non-Gaussian wave statistics are comparable to 
those existing in the second order of scattering from a normally distributed sea. A better understand- 
ing of the measurement problem - to this level of detail - will clearly require muih additional analy- 
sis work and experimentation. 

In a companion paper (Jackson et al. t 1980, in this issue), experimental aircraft data are pre- 
sented that support the basic findings of the present theoretical investigation: namely, that direc- 
tional wave slope spectra can be measured with rather good fidelity from near-vertical measurements 
of microwave backscatter, remotely, using a scanning beam, short-pulse radar system. 
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APPENDIX A: Analysis of SP and DF Systems 
Consider the SP system of Fig. 3a. If the FT of the backscattered field is given by Eg(v) ■ 
then the FT of the square-law detected signal is 

/* 0 ® 

P(w)- S(*)S*(i'-t.>)E o (i0E$(v-w)dv (Al) 


The ensemble average output of the spectrum analyzer is the pulse spectrum P(w) ■ <|/*(co)| 2 > as 
seen through the Hj bandpass filter; i.e., P out (< 0 |) ■ (2 a) 3 / |H ( (w)| 2 P(u>)d<o. If the analysis band- 
width 0 ti of the Hj filter is small compared to the variation in P(w), then P out (wj) “ 0 ai P(U|). From 
( Al ) we have 


P(u>) 



M(v, v\ w) E 0 (i») E* (a - w) E* («>') E 0 (v' - u)dvdv' 


(A2) 


where 


M = <S(v) S*(v - w)S*(v’)S(v’ - w)> 


(A3) 


The moment M is given by (55) and (56) 

M-(V f 57/T r )t<|S(v)| 2 >l 2 * |5(w) + 6(Av) + P mod (w) + P inod (A»-)} (A4) 

where Av ■ v - v\ For weak modulation. /P mod ( Ai^)dAv « 1 , and we can neglect the contribution 
of P mod (Av) to P(cj); since (|S(v)| 2 ) is practically a constant over the few percent pulse bandwidths 
we are considering, we then have on substituting (A4) into (A2): 

P(w)*(>/J7/T r )l<|S(v)|2>l 2 {l/E o (*0E>-cu)d*| 2 • 

• l*(w) + P rood (w)l ♦/ E 2 (i»)E 2 (v- w)dv} (A5) 

Pulse waveforms are often Gaussian in shape; iet the FT of the envelope of the incident waveform 
be given by 

E 0 (v) » exp£- -i- (*/flp) J 


(A6) 
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Then P(w) has the form 


P(w) * (\/H/T r ) 1<|S(*0I 2 >] 2 



| 6(W) + y/2irfi p + P|nod (w, ^| 

Now, consider the DF system diagrammed in Fig. 3b. The and Hj are non -overlap ping bandpass 
fflters with center frequencies » i and v^, respectively; the K t and Kj are dc blocking filters. The out- 
put of the correlator can be written as 


P out ("i-V* 2 *) 3 J K i (u > )K;( W )<P i /7>d W 

(2«) 4 JJ M(i», v\ u) E 0 (v) • 

• -w)Hl > (v’) E*(«»')H j (/-w)E 0 (»/'-w)d^di'' <A8) 

where M is given by (A4). Since the K filters block dc, Kj(0) ■ Kj(0) * 0; let Kj « Kj ■ (2er 1 , for 
w> 0. Then, assuming that the H-filters are non-overlapping and that the bandwidths 0 gi and 0 t j 
are narrow compared to the variation in P m(N |, we have using (A4), approximately, 

P ou* * (2ff ) 5 P mod^i " V * J HtfAu (A9) 

where //j(w) * J Hj(t') E o (t0 w) E*(i»- toldHi -♦j). Two waveforms of interest are SP wave- 

forms and DF waveforms that are matched to the H-filters. For SP waveforms. E 0 is constant over 
the bandpasses of the Hj filters. If the filter functions are Gaussian with standard deviations 0 ai and 
0 aj . then 

P ou«“0.u P <*V (A,0) 

where and 0 aij * 0 ai 0 a j (20^ + 20 a j)"^. A similar analysis may be carried out for the 

case where the receiver filters are matched to the transmitter waveform, e.g., if the DF waveform is 
t 0 * t oi + E oj then H, » (i - j). 


APPENDIX B: Alternative Derivation of PO Limit 


The fact that Eq. (7) has the unique stationary points given by (8), of course, has to do with 
the smoothness and extremal properties of the mathematical expectation, i.e., the factor in braces. 
The following development may clarify this point Consider the second moment of S(k) in the case 
of backscatter in two dimensions only. Then, we must evaluate the following integral in the limit 

k 4 

• f 

I fc 3 2kcos0 r < e - i2kc ose ( r(w)-r(0)l >e -i2ksinew dw (Bl) 

«/dOO 

Let z = [ £(w) - £(0)1 /w and let p z (z; w) be the pdf of z. Then 

-i» 

4 

PO 

I k = 2kcos0 JJ p z (z; w) e i2kcosa(z-taii 9 )w dzdw (B2) 

.OO 

Since p z is a slowly varying function of z and w, we can apply the method of stationary phase. A 
saddle point exists at z = tan0 and w = 0. Letting 

u = (z - tan0 + 2kcos0w)/\/T 
v = (z - tan0 - 2kcos0w)/\/T 


then since limz = £ where £ is the surface slope and since lim p z = p t it follows that 


w-»0 


w-*0 


we 

lim I k = p fx (tan0) • I* f e ‘l(u 2 -v 2 )/2] dudv 
f-frOO J J 


k-*o° 

= 2jt p rx (tan0) 

which result agrees with Barrick ( 1 968b). 


(B3) 
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APPENDIX C: Socondordor Terms in the Expansion of M[(K) 
The second-order terms in the expansion of ^ , eq. (38), are: 

- | [(Kcotei 4 R : + R. ( R „ ■u.u, v„v, + 

+ (Kcotfl) 2 R a R* (u tt -2u o v () + v 0 v fl ) - 
- 2(Kcot0) 2 R R >a<J u 0 v fl + 2(Kcot0) 3 R R 0 (u a - v # ) - 


(Cl) 


- 2(Kcot0) R a<| R y u a y fi (u y ~v,)j 

Fourier transforming with respect to u and v according to eq. (21 ) and making use of the identities 
(3d), we get for the second-order terms in the expansion of S, cq. (40): 

-l-(Kcotfl) 4 R 2 + ~ R_.R 


J 

P 2 


,«0 ,*« 


R,. R,,-(Kco<«) ! ^ P - R R„ 
* :i(KootO) 5 R R,„ - iRcotO R „ R , 


(C2) 


f rom the definition of the height spectrum as the Fourier transform of the covariance function, vi:.. 


(K) a FT{r} a — R(w)e” i ~* - dw 

Cff) 2 J 


(C3) 


it follows that the FT of the second-ordsrBertm in S is given by various convolutions involving the 
height spectrum; e.g. ( 
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FT{R R J - FT{r} • FT{R >U } 

* F(K) • (iK a F(K)1 

* J°° F(K') [i(K; - K tt ) F(K* - K)1 dK’ 


(C4) 


If an overbar is used to distinguish a variable wavenumber entering into a convolution, then the 
second-order terms in the expansion of P mod , eq. (48), are written as 


4 (Kcotfl) 4 F • F + i K.K.F * KK. F 

2 2 n 2 * y t 


1 P,*« P,0y 


, / p,«(» p,« p, A — — 

+ (KcotO) 2 K a F*K fl 


F + 


+ (KcotO) 2 


P,a Pj 


F * K K. F - 2(Kcotd) 3 — F * K a F 


,2 • 0 ’ p 


- Kcot0 


^ P,a7 P.g*P.gy P,^ 


K tt K fl F*K t F 


(C5) 
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FIGURE CAPTIONS 

Fig. 1 . Overall geometry and coordinate definition. 

Fig. 2. Simple tilt model of reflectivity modulation. 

Fig. 3. (a) Schematic of modulation noncoherent short pulse system, (b) Diagram of wideband 

dual frequency system (after Parzen and Shiren, 1 956) 

Fig. 4. Possible realization of a modulation coherent short pulse system, r denotes fast signal 
time and t slow signal time. The slow time filter can be realized by a sample-and-hold 
that is triggered at successive delays given by ft, followed by an accumulator. The spec- 
trum analyzer can be realized by digitally fast Fourier transforming the slow time 
filtered pulse return, and squaring the magnitude. 

Fig. 5. Sample calculations of the spectrum P mo( j(K) of backscattered impulses from a Gaussian 
sea surface in two dimensions in the second order of scattering, (a, b) The surface slope 
spectrum K 2 F(K) is the Phillips spectral form, BK" 1 , with an abrupt cutoff at K 3 K,,. 
The slope spectrum and P mod are both normalized by their values at K. c . The calculations 
are performed for a range of wave steepness 5 Q = K Q o 3 n/B/ 2, incidence angle 0 and wind 
speed U (which determines the total mean square slope according to eq. (30)). In all cal- 
culations, the location of the peak at K Q is preserved; low frequency whitening - or inter- 
modulation (IM) - is exhibited; and the frequency response (FR) at wavenumbers K> K Q 
is in the form of a droop. Panels (c-e) show the variation of the IM at K 3 0.5 K 0 and the 
FR at K = 2 K 0 over the range of parameter values indicated. 

Fig. 6. Spectrum of backscattered power when a short pulse of bandwidth 0 p is transmitted, and 
definition of the SNR. 

Fig. 7. Post-detection ‘‘pulse-Doppler” spectrum. 
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Table I. Parameter Values Used in Sample Calculations 


Incidence angle 

0® 12.5° 

Slant range 

r 0 ® 700 km 

Satellite velocity 

V= 7 kms" 1 

Beam spot size* in range 

L « 5.0 km 

p 

Beam spot size* in azimuth 

L. - 2.5 km 

EM wavenumber 

k® 2»/2cm 

Water wrvenumber 

K ® 2 j /200 m 

Wind speed 

U® 18 ms -1 

Pulse bandwidth parameter^ 

0 p = 1.4 X 10 8 

Surface range resolution parameter 

Ap = 5 m 


*One-*igma value of Gaussian, one-way power pattern 
t One-sigma value of Gaussian envelope of field strength 
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X! (REFERENCE) 


Q 


surface range 


The fractional backacatter cross-section variation is 

8o_ 6 ( q ° A) , W 
0 a°A " o° + A 


Assume that o° ■ i r sec 4 © p(s) where s ■ tan© toc and p is the slope pdf. 
From the spherical triangle, if ©»0, it follows that 0| OC ~ © - 0 cosa; 
or ©| OC ~ 8 - p*vf. Since A « cr/2sin© )0C it follows that 
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